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(57) Abstract: The present invention relates to classification and lo- 
calization of arrhythmia* More specifically, a system and method is 
provided for oevetoping a database of boo^ surface ECO P wave maps 
for classification and localization of left-sided atrial arrhythmia. The 
invention includes generating and receiving P wave data in a subject 
by left atrial paring or recervingP wave da^ 

nmeousiy occurring or induced left atrial arrhythmias; computing (e.g. 
potential or mtegraOniar« of the P wave o^^ 
specific to a left atrial ectopk origin; verifying the classification pro- 
cedure; averaging the classified maps into mean maps; and storing and 
accessmg the mean raai>s in the database. The in^ 
data in the database can be used to automatically classify and localize P 
wave data from a patient obtained c^ng a left atrial anliythm^ 
atrial tachycardia, focal atrial fibrillation or orthodromic atrioventricu- 
lar reentrant tachycardia. 
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HELD OF THE INVENTION 
This invention relates generally to a method for classification and localization of 
arrhythmias. More particularly, the present invention relates to a system for 
developing a database of body surface electrocardiogram (ECG) P wave data maps for 
tecafizadon "f «he origin of left-aided atrial arrhythmias. 

BACKGROUND ART 
Early investigations utilizing atrial pacing in man to stndy the variations in scalar ECG 

morphology and vector loop of the P wave with varying sites of focal atrial origin, 

were primarily geared towards die differentiation of left- from right-sided atrial 

rhythms {e.g. Massumi RJL and Tawakkol AA. (1967), Direct study of left atrial P 

waves, Am. J. Cardiol 20:331-340; and Harris B.C. et at (1968) Left atrial rhythm; 

Experimental production in man. Circulation 37:1000-1014). Despite several attempts 

to develop a set of morphologic ECO criteria specific to a left-sided origin of an 

ectopic atrial rhythm, different algorithms were proposed and consensus appeared 

difficult to attain. 

Mirowski in 1967 initially suggested that a negative P wave polarity in lead I and a 
"dart and dome" configuration in V, was related to a left-sided origin and later added 
that a negative P Wave in V 5 was more specific particularly when the aforementioned 
features were absent {Mirowski ML (1967) Ectopic rhythms originating anteriorly in 
the left atrium; Analysis of 12 cases with P wave inversion in all precordial leads. Am. 
Heart J. 74:299-308). Although these findings were partly underlined by others (e.g. 
Leon D.F. etal (1970) Right atrial ectopic rhythms; Experimental production in man. 
Am. J. Cardiol 25:6-10), Harris et al. in 1968 (same reference as mentioned above) 
contested the importance of P wave inversion in leads I and V 6 for a left-sided rhythm 
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and alternatively proposed that a terminally positive P wave in V, was a more specific 
finding for a left atrial origin. Conversely, Massumi and Tawakkol (same reference as 
mentioned above) noted a profound variability in P wave morphology after electrically 
sttrmilating the left atrium from similar anatomical sites in different patients and did 
not believe that distinct BCG criteria specific to areas of ec^ 
formation could be developed. 

Usmg temporarily implanted pacing wires following cardiac surgery, Maclean et aL 
sobsequendy perform^ 

at a total of 12 epicanfial regions {Modem WA.H. et al (1975) P waves during, 
ectopic atrial rhythms in man; A study utilizing atrial pacing with fixed electrodes. 
Circulation 52:426434). Overall the results of this study were disappointing in that 
only a few specific correlations between P wave morphology and site of origin could 
be made: (1) a negative P wave in the inferior leads with pacing of the inferior regions 
in either atrium; (2) a negative P wave in lead I with left atrial pacing near the left 
pulmonary veins; (3) a positive bifid P wave in V, with pacing near the lower 
pulmonary veins and coronary sinus. Therefore, these investigators and recently 
others using contemporary mnltisite endocardial catheter mapping techniques (Man 
K.C etal (1996) Spatial resolution of atrial pace mapping as determined by uniploar 
atrial pacing at adjacent sites. Circulation 94:1357-1363), underline the complexity of 
visual analysis of the low-voltage P wave and concluded that the 12-lead ECG was of 
limited clinical value in localizing ectopic atrial foci. 

The hmited role of the standard 12-lead BOO to distinguish right from left atrial focal 
activity has also been reported by Kahnan et aL during ablation of atrial tachycardia 
guided by intracardiac echocardiography (Kalman J.M. et al (1998) -Crista! 



3 



WO01/7«« FCT/USDl/mW 



tachycardias": origin of right atrial tachycardias from the crista terminate identified 
by intracardiac echocardiography, J. Am. Coll. Cardial 31:451-459). Kalman et al. 
reported that the P wave morphology on the 12-lead ECG of tachycardias arising from 

the upper part of the crista temrinalis in the ri ght atri m n or the right upper pohnonaty 

ft*/ £■< - . - *..« y/y'* («n9- : ? ; -• ■ 

vemm the left atrium, demcastr^ 

structures. TangetaL suggested previously that achange in P wave morphology in 
lead V, from Diphasic during sinus rhythm to completely positive during tachycardia 
wouM be helpM in discriminating foci arising from the superior portion of the crista 
terminate and the right upper pulmonary vein {Tang CW1 et al (1995), Use of P wave 
configuration dming atrial tachycard^ 

26:1315-1324). The latter finding was found to be predictive for right upper 
pulmonary vein tachycardia foci even when lead aVL would demonstrate a positive 
instead of a negative P wave. 

There have also been attempts to localize the atrial insertion site of an accessory 
pathway during orthodromic atrioventricular (AV) reentrant tachycardia. Farshidi et 
at initially demonstrated that a negative retrograde P wave in lead I was associated 
with a left-sided atrial insertion of the accessory pathway (Farshidi A.et aL (1978) 
Electrophysiologic characteristics of concealed bypass tracts: clinical and 
electrocardiographic correlates, Am. J. Cardial 41:1052-1060). Garcia Rivera et al. 
was able to separate free wall accessory pathway locations in the left and right atrium 
using the retrograde P wave polarity in leads I and V, (Garcia Rivera Ketal (1980) 
Retrograde P wave polarity in reciprocating tachycardia utilizing lateral bypass 
tracts, Eur. Heart J. 1:137-145). Other investigators have later studied the localization 
resolution of the 12-lead ECO for this particular application by confining pace 
mapping to the annular regions of the left and right atrium (e.g. Fitzgerald DM. etal 
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(1996) P wave morphology during atrial pacing along tJxe atrioventricular ring; ECG 
localization of the site of origin of retrograde atrial activation, J. ElectrocardioL 29:1- 
JO). However, despite a more directed pace mapping approach, the paced P wave 
morphology only allowed a gross separation of pacing sites in terms of a left- versos 
right-sided origny an inferior origin m either atrium, Aran nrigin fn the right fee walL 
In contrast, Tai et aL recently reported that the polarity of the retrograde P wave in 
leads I,n T m, aVF, and V x obtained during AV reentrant tachycardia allows accessory 
pathway localization to 9 possible annular regions with an overall accuracy of 88% 
provided a clearly visible P wave could be dtscriminated gbf CH etjaL (1997) A new 
electrocardiographic algorithm using retrograde P waves for differentiating 
atrioventricular node reentrant tachycardia from atrioventricular reciprocating 
tachycardia mediated by concealed accessory pathway, J. Anu ColL Cardiol 29:394- 
402). It has also been realized with this application of the surface ECG, that visual 
analysis of the low-voltage retrograde P wave is frequently hampered by the 
simultaneous occurrence of die preceding cardiac cycle's high-voltage TU wave. 

As mentioned above, visual inspection of the 12-lead ECG in an attempt to correlate 
changes in P wave morphology and polarity with various locations of ectopic activity 
is hampered by a low-amplitude signal which is often buried in the preceding TU 
wave, and also by inter-patient inconsistency in P wave pattern with comparable 
locations of ectopic left atrial activation, and poor overall spatial resolution in both die 
left and right atrium. It has been recently demonstrated using right atrial pace mapping 
that a significantly higher electrocardiographic localization resolution is clinically 
feasible when multi-lead ECG techniques are adopted (SippensGroenewegen A. etaL 
(1998) Body surface mapping during pacing at multiple sites in the human atrium; P 
wave morphology of ectopic right atrial activation, Circulation 97:369-380). 
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There is one preliminary clinical report investigating the use of multi-lead ECG or 
body surface mapping to discriminate ectopic left atrial activation (Kawano S. and 
Hiraoka MJ>- (1995) P wave mapping in ectopic atrial rhythm, bu Yasui 5. et aL Eds. 
"Advances in Body Surface Mapping and High Resolution ECG": Nagcya, Japan: 
Ufe Medacom pp. 47-56). Kawano and ffiraoka showed that endocardial pacing at one 
right 0ow) and three left flow, middle, and high) atrial locations produced 
characteristic P wave body surface potential map patterns. Access to the latter left 
atrial locations was attaint However, only a 

moderate number of pacing site locations were studied while pacing from the 
p ulmo nary veins was not attempted. 

It is clear thai the development of noninvasive methods for clinical localization of left 

atrial arrhythmias with sufficient and high electrocardiographic resolution is difficult 

Therefore, noninvasive localization of left atrial arrhythmias is in need of alternative 

methods that enable improved discrimination and differentiation of 

electrocardiographic P wave patterns or morphologies. 

OBJECTIVES AND ADVANTAGES 
In light of die above* it is the primary objective of the present invention to provide a 

system for developing a database of body surface ECO P wave data maps. Body 

surface ECO mapping can be applied toimprove noninvasive localization of left atrial 

arrhythmias prior to radiofrequency catheter ablation. 

His an additional object of the invention to provide a method of storing a database of 
mean body surface ECO P wave data maps in a computer-readable medium for 
retrieval, analysis and sharing. 
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11 is an additional object of the invention to provide a database that will enable detailed 
noninvasive localization of left atrial arrhythmias so that invasive catheter mapping 
may be targeted and minirmzed to a sinaflcoiulnedseginent of the left atrium prior to 



The advantage of the present invention is that the present invention resn^ 

and more efficient mapping and ablation procedures with less fluoroscopic exposure of 

the patient and physician. Additionally, the patient will undergo less trauma and 

bene&fomiredoction m the risk for dot fonnatira 

for catheter manipulation in the left atrinmcanbedgnfficantiyieduced. 

An additional advantage is that the database can be used in an outpatient clinic or 
physician's office to noninvasive^ screen patients and map the arrhythmia's origin 
before admitting a patient to the electrophysiology laboratory to undergo ablation 
therapy. This noninvasive diagnosis and screening will also enable improved decision 
making regarding the optimal choice of therapy, e.g. anti-arrhythmic drags, permanent 
pacing, ablation treatment, or combinations of these treatment modalities. Another 
advantage of the present invention is that it wffl also allow for noninvasive follow-up 
of patients with left atrial armythnrias after therapy delivery to assess arrhythmia 
recurrences due to treatment failure or disease progression. 



SUMMARY 

The present invention provides a method for classification and localization of 
aniiythmias. More specifically, the present invention provides for a system for 
developing a database of body surface ECG P wave data maps for classification and 
localization of left-sided atrial arrhythmias. 
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In accordance with exemplary embodiments of the present invention the method 
generally comprises receiving P wave data from two or more subjects; computing 
maps of the P wave data; classifying the maps and thereby creating classified maps; 
verifying the classified maps; aveiagmg the classified and verified in^ 
cnatmgrneanmarg; and storing and m^^g the mean naaps in the database. 

The P wave data may generally be obtained in a subject by electrically s timu la ting or 
pacing the left atrium using a probe by fiw example a retrograde aortic or a rranssep^ 
procedure. Alternatively, the P wave data may also be obtained during a 
spontaneously occurring or induced left atrial arrhythmia. 

The maps may be computed from the P wave data in the form of P wave potential 
maps or P wave integral maps. The maps are classified based on spatial pattern 
characteristics. 

The maps may be verified before averaging. The mean maps may then be calculated 
from the verified maps to construct a database of characteristic body surface ECG 
mean P wave maps specific to the origin of left atrial pacing or left atrial arrhythmia. 
The mean maps of the P wave data in the database can be used to automatically 
classify and localize P wave data from a patient obtained from a left atrial arrhythmia. 
For instance, P wave data can be obtained from an atrial tachycardia whereby the 
origin is localized, a focal atrial fibrillation whereby the origin of a single atrial 
premature beat or an initiating beat is localized, or an orthodromic AV reentrant 
tachycardia involving an accessory pathway whereby the atrial insertion of an 
accessory pathway is localized. 
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BRIEF DESCRIPTION OF THE FIGURES 
The objectives and advantages of the present invention will be understood by reading 
ttefolkw^ de ta iled dcsc x qitioa in conjunction with the drawings, in which: 
PIG. 1 depicts an overview of the development of a database of mean body surface 
ECG P wave data maps for classification and localization of left atrial 
anhyt h mias according to an embodiment of the present invention. 
FIG> 2 depicts a system for developing a database of mean body surface ECG P 
wave data maps for classification and localization of left atrial arrhythmias to 
better classify and localize left atrial arrhythmias in patients according to an 
embodiment of the present invention. 
FIG. 3 depicts a method for separating obscured P waves from a superimposed 

QRST complex according to an embodiment of the present invention. 
FIG. 4 depicts a method for selecting and processing good quality signals according 

to an embodiment of the present invention. 
FIG. 5A depicts a sensor system having an array of sensing locations distributed 
across a subject's torso surface according to an embodiment of the present 
invention. 

FIG. SB depicts a body surface P wave integral map; a plot of a data matrix 
generated by mapping the integral values with positions corresponding to the 
location of the sensors across a subject's torso surface according to an 
embodiment of the present invention. 

Fig. 5C depicts an ECG tracing illustrating a method for calculating an integral 
value across a selected rime m 

location according to an embodiment of the present invention. 
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Fig. 6 depicts six mean body surface P wave integral maps obtained during pacing 
at the left upper or left lower pulmonary vein according to an embodiment of 
die present invention. 

KG.. 7 depicts a database of 17 different mean P wave integral maps generated by 

left atrial pacing according to an embodiment of die present mveniion- 
FIG. 7A depicts different mean P wave iQtegndniaps and n^ associated segnients 

of pacing origin represented in a postero-anterior view of the left atrium 

according to an embodiment of the present invention. 
FIG. 7B depicts different mean P wave inte 

of pacing origin represented in an antero-posterior view of the left atrium 

according to an embodiment of the present invention. 
FIG. 7C depicts different mean P wave integral maps and their associated segments 

of pacing origin represented in a left posterior oblique view of the left atrium 

according to an embodiment of the present invention. 
FIG. 8 depicts the design of an apparatus for classifying and localizing left-sided 

atrial arrhythmias according to an embodiment of the present invention. 
FIG. 9 depicts a block diagram of storing a database of mean body surface ECG P 

wave data maps in a computer-readable medium for retrieval, analysis and 

sharing according to an embodiment of the present invention. 



DETAILED DESCRIPTION 
Although the following detailed description contains many specifics for the purposes 
of illustration, anyone of ordinary skill in the art will readily appreciate that many 
variations and alterations to the following exemplary details are within the scope of the 
invention. Accordingly, the following preferred embodiment of the invention is set 
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forth without any loss of generality to, and without imposing limitations upon, the 
claimed invention. 

The present invention involves a system for developing a database of body surface 
ECOP wave data maps which allows a more mtuitrve and direct m^ 
heart signals in terms of the under^caolac source. This database enables detailed 
noninvasive localization of left atrial arrhythmias by correlating subtle spatial 
differences in P wave morphology with various locations of ectopic atrial activity. 
71* database is derived ftom P w 

development be updated using additional and/or future P wave data of subjects to 
improve its resolution for discriminating focal atrial activity. 

Clinical application of the database of P wave data maps enables improved 
noninvasive localization of left atrial arrhythmias prior to radiofrequency catheter 
ablation. Noninvasive localization of left atrial arrhythmias is important since it can 
target and nrinimize invasive catheter nnipping to a small confined segment of the left 
atrium prior to ablation therapy delivery. Consequently, shorter and more efficient 
mapping and ablation procedures wiD result in less fluoroscopic exposure of the 
patient and physician as well as the patient will undergo less trauma and benefit from 
reduction in the procedure-related risk for clot formation and stroke. The database of 
the present invention can also be applied to noninvasively screen patients and localize 
armythmias before admitting a patient to the electrophysiology laboratory to undergo 
ablation therapy. This noninvasive screening and arAythmia localization can easily be 
performed in an outpatient clinic or physician's office and enables improved decision 
making regarding the optimal choke of therapy, eg. anti-arrhythmic drugs, permanent 
pacing, or ablation treatment, or combinations of these treatment modalities. 
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Additionally, the database of the present invention also enables noninvasive follow-up 
of a patient to assess arrhythmia recurrences due to failure or disease progression after 
therapy delivery has taken place. 

The embodiments of this invention relate to developing a database of mean multi-lead 
body snrface P wave data map patterns mat can be used to classify and localize left 
atrial arrhythmias. Examples of currently known left atrial arrhythmias are atrial 
tachycardia, focal atrial fibrillation, and orthodromic AV reentrant tachycardia 
involving an accessory pathway. The present mverrtion may also be applied to other 



The use of body surface mapping allows ECO data presentation in a unique spatial 
format that contains a display of the P wave distribution on the entire torso surface. 
This electrocardiographic display format provides a better way to identify the focal 
source or reentrant circuit of atrial tachycardia, the focal source of atrial fibrillation, as 
well as the atrial insertion site of an accessory pathway with orthodromic AV reentrant 
tachycardia. 

The method of the present invention is shown in the simplified flow diagram of FIG. 

Ic* ■•,*•>-• i5 lr 

1 The electrical signal from a heart cycle 2 is a spontaneously occurring or induced 

s\ •'V-v *" 

left atrial arrhythmia or it can be obtained using left atrial pacing. One or more heart 
cycles 2 can "be obtained in two or more subjects: The heart signal is received 4 as 
unipolar or bipolar ECO signals simultaneously obtained from an array of multiple 
electrodes on piedeteiriuned torso sites. During signal processing 6, leads with 
inferior signal quality are rejected and baseline adjustment is performed by linear 
interpolation. A series of potential maps of the P wave cycle is computed at 2-ms 
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intervals (other intervals can also be used). In a potential map, the signal sensed at a 
time instant by sensors placed on a patient's torso is plotted on a renngmtefogof ; ite 
torso, e.g. as a contour map or color map. P wave onset and offset are determined, for 
L instaiK*, at the tinie instant when a p^ 
when tte peak vrftage drops bi^ Other criteria \ v 

as well as an automated detection of P wave onset and offset can also be used. ' ^ 
SubsequenUy, a body surface ECO P wave integral map is computed 8 over this 
selected time interval Ih an integral map, the signal seijsed by sensom is mtegrated 
overasriectedthasmtervalandt^ 

e.g.asaamtcmrn^orcolormap. QoaKtative data classification 9 is carried oat by 
grouping m aps based on pattern similari ty. Other methods of data classification based 
on, e.g. pattern recognition algorithms, neural networks, statistical routines or 
anatomical location of the left atrial ectopic origin, can also be used. Mathematical 
analysis using eolation coefficients is then performed to verify the qualitative data 
classification procedure. FurthemK>re,j3ean^a^^ D f the classified 

maps. The mean maps (also called mean P wave maps) obtained during left atrial 
pacing or spontaneous or induced left atrial arrhythmias are stored in ajatabase 12. y f f 
this unique databa se will then conta mmean body surface P wave^ta'maps for 

andJocaBzation of left atrial arrhythmias.. The n^ i 
reference for localizing P wave data from a patient with left atrial arrhythmias such as 
atrial tachycardia, focal atrial fibrillation, or orthodromic AV reentrant tachycardia % ' 
using an accessory pathway. By ccrriparing a patient's arrhythmia data to the maps in 
the database, a physician can provide a better diagnosis and make better decisions for a 
treatment regimen for an individual patient. 

\J y 13 •:• 
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Referring to FIG. 2, there is shown by means of a more detailed general flow diagram 
die overall sequence of steps to develop a database of body surface ECGP wave data 
maps for classification and localization of left atrial arrhythmias. One or more heart 
cycle sigpafc 20 am acquired in tro The electrical heart cycle signal 

20 is obtained during left atrial pacing or spontaneously occurring left atrial 
arrhythmias Left atrial pacing is performed using an electric stimulus probe such as a 
catheter using a transseptal or a retrograde aortic procedure. Alternative approaches or 
techniques other than a catheter are also feasible. Spontaneous left atrial arrhythmias 
can fcrexan^Je be atrial premature beats, atrial tachycardia, to 
atrial fibrillation, and die retrograde atrial activation during orthodromic AV reentrant 
tachycardia involving an accessory pathway. Left atrial arrhythmias can also be 
induced by programmed stimulation, pharmacologic challenge, or other means. 

■"HjdivJta 

In a particular embodiment, left atrial pacing is conducted at 79 pacing sites (higher \j.r ' 
and lower numbers can also be used) under biplane fluoroscopic guidance in the 45 ^ 
degree right and 45 degree left anterior oblique projections. Other fluoroscopic 
projections or other imaging techniques may also be used. A pacing protocol is aimed 
at accumulating pacing sites in key anatomical structures (Le. pulmonary veins, around 
the circumference of die mitral axmulus through die coronary sinus or directly from 
within the left atrium, the left atrial appendage, and the septum) as weU as at obtaining 
additional coverage of the remainder of die left atrium. All pacing positions on die 
biplane fluoroscopic images are visually translated into positions cm an anatomical 
diagram of the left atrium while die pacing protocol is in progress or immediately after 
its completion. 



14 



WO 01/76461 



PCT/US01/12263 



The heart cycle signal 20 is processed by a P wave classification routine 21. At step 
22, the heart signal is received as unipolar or bipolar ECO signals simultaneously 
obtained from an electrode array on predetermined torso sites. Referring to FIG. 5A, 
the techniques of the present invention will generally make use of an array 100 of 
sensors 112 distributed across anterior and posterior skin surfaces of torso T on patient 
P. Array 100 provides mum-lead ECO data at a plurality of sensing locations 
distributed across the torso T, typically at over 20 sensing locations, more preferably at 
over 40 sensing locations, and ideally at 62 or more sensing locations. 

Sensors 112 generally comprise mupolar or bipolar electrotfcs coupled to the patient's 
skin, or to an alternative accessible body surface (for example via a transesophageal 
approach) suitable for measuring electrical body surface potential. Suitable electrode 
structures may include those described in U.S. Patent Nos. 5,311,873 and 5,634,469, 
which are incorporated herein by reference. Alternative sensor array structures and 
associated data acquisition and manipulation components may also be used. 

ECO data is preferably acquired simultaneously from each of sensors 112 of array 100 
at a sampling rate of over about 500 Hz, ideally at a sampling rate of about 1,000 Hz 
or more. In some embodiments, sequential sampling of sensors 112 from array 100 
may alternatively be used, and higher or lower sampling rates also are feasible. When 
a lower sampling rate is used, the data may be upsampled using multirate filter banks. 

Referring to FIG. 2, at step 24, if the P wave is obscured, e.g.. due to the P wave 
activity coinciding with ventricular recovery of the preceding cardiac cycle, the QRST 
subtraction 26 routine can be invoked to separate the signals. It may be possible in 
some circumstances to artificially separate die P wave from the QRST interval, i.e., by 
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the use of pharmacological compounds that temporarily block AV conduction or by 
electrical pacing of the ventricle. Alternatively, other signal separation methods and 
systems may also be used. 

Refemng to FIG. 3, the QRST subtract 

constructed from averaged QRST complexes obtained from body surface ECG 
measurements to enable isolation of toe otherwise obscured P wave activity. 

The method of 26 generally includes itcordiiigof un^olarorb^larBCGdatefiOTn 

me array of torso sites in steps 50 and 54, The signal 50 is obtained during left atrial 

pacing or during a spontaneous or induced left atrial arrhythmia that contains electrical 

heart signals in which atrial activity is obscured by ventricular activity. This recorded 

signal 50 will also include both the P wave and a superimposed QRST signal portion 

52 and is further referred to as the obscured P wave signal. The method 26 also , . 

receives at least two template signals 54that contain electrical heart signals obtained ^ 0 /^ /a j. 

during normal sinus rhythm or atrial overdrive pacing for the construction 56 of a 

QRST template 57. fa die exemplary embodiment, about 100 cardiac cycles of 62- 

channel unipolar ECO data are recorded, although fewer cycles may be used if the 

spatial and temporal variations of the QRST complex are relatively low. Typically, 

more man ten cycles will be used, often more than 50 cycles for construction 56 of the 

QRST template 57. 

An adaptation of the QT interval of the QRST template can take place to correct for 
differences in heart rate with the obscured P wave signal. One or mare fiducial points 
and windows are identified and annotated 58 for boto the QRST teinplate to 
called annotated QRST template as well as for the obscured P wave signal, to create a 
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so called annotated obscured P wave signal. Alignment 60 of the one or more fiducial 
points and windows annotated in the annotated QRST template and the annotated 
obscured P wave signal then takes place. After the alignment 60 of the annotated 
QRST template and ann ot ate d obscured P wave signal, the annotated QRST template 
is resampled and modulated 62 to fmte compensate for remaining discrepancies in 
du ra tion and/or voltage. This creates a resampled and modulated QRST template. The 
atrial activity contained in the obscured P wave signal is isolated by subtracting 64 the 
resampled and modulated QRST template from the obscured P wave and thereby 
generating an isolated P wave 66. It is important to realize that the QRST template is 
adaptive to optimize the subtraction by first aligning the QRST template and then 
resampling and modulating the QRST template. The reason is to account for a 
discrepancy in heart rate and voltage amplitude of the QRST template with the TU 
complex of the obscured P wave signal. The method 26 includes a variety of 
techniques that can be used to modulate and resample the QRST template. Generally, 
this approach of modulation and resampling allows for the surface EGG measurements 
to retain their intricate spatial and temporal detail within the P wave morphology. The 
method 26 is capable of unmasking and preserving subtle electrical heart signal details 
of the relatively low-voltage P wave despite the obscuring superimposed relatively 
high-voltage QRST complex. 

Further details of QRST subtraction method 26 are described in. and are incorporated 
herein by reference, U.S. Patent Applications Serial Nos. 09/809,719 and 09/808,728, 
both filed March 14, 2001. 

Referring to FIG. 2, all lead tracings are screened at step 28 in order to reject poor 
quality signals or signals from read sites that are dislodged or obscured by 
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defibrillation patches on the front and back of the subject's chest Hie details of the 
agnal-piocessing step 30 are depicted in the flow diagram of FIG. 4. All lead tracings 
are visually or automatically screened for signal quality 70. Absent signals 72 are 
detected and deleted 74, and substitute values 76 are interpolated based on adjoining 
lead recordings. Poor quality signals 77 are detected and rejected 78» and solwtitote 
values 80 are mterpolated based on adjoin^ mterelectrode off set 

differences 82 are corrected 84 and Jiiiearbaselme drift 845 is corrected 88. Both the 
offset and drift may be corrected 84, 88 by using linear interpolation after selecting 
one or more isoelectric time instants or intervals. Digital filtering to remove line 
frequency interference 90 is performed when necessary with a filter 92 such as a 60- 
Hz notch filter or other type of filter. Other means of filtering can also be used. The 
end result of the signal processing 30 is a good quality signal 94. 

Referring to FIG. 2, by way of example, P wave integral maps are computed 32 as 
follows, mectrocardiographic data is analyzed by determining the onset and offset of 
the P wave. Subsequently, an integral map is computed over this selected time interval. 
Referring to FIG. 5C, graph 114 includes an ECG signal tracing 116 representing the 
variation in voltage over time. Signal tracing 116 may be used to evaluate heart cycle 
signals. In general, one or more reference heart cycles will be selected for 
manipulation and comparison. ECXJjracjng 116 can be used to fetermine an onset 
(starting time) 118 and offset (ending time) 120 of a time interval 122 of the heart ^ \ .. 
signal cycle that is of interest for evaluating one or more regions of tire heart. For 
instance, P wave onset and offset can be determined at the time instant when a peak 
reference voltage, e.g., 30 jrV is reached and when the peak voltage drops below the 
reference voltage, respectively. Alternative criteria for determining P wave onset and 
offset might also be used, and automated drtection of time interval 122 is also feasible. 
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Instead of a time interval, single or multiple time instants in the P wave cycle may also 
be selected 

FIG. 5A shows an exemplary map of sensor locations 112 on the front and back of a 
subject's torso. FIG . 5B shows an exemplary plot of a data matrix generated by 
mapping the integral values with positions conespomlin^ 

112 across a torso T of a patient Pin FIG* 5A. FIG. 5C shows the calculation of an 
integral value from a plot 114 of heart cycle signal 116. In the exemplary embodiment 
of FIG- 5G, a P wave integral nu me ri cal value may he calculated based on heart cycle 
signals 116 within selected time interval 122 having an onset 1 18 and an offset 120 for 
a particular sensor location N10. This calculated P wave integral value reflects the 
time/amplitude area of ECG signal at the sensor location within the selected time 
interval. Similar integral values are calculated for each sensor location, and the sensor 
values are mapped within data matrix 124 continuously from a portion of the data 
matrix associated with a front F of torso T, across a side of the patient P, and to a back 
B portion of torso T\ As shown in FIG* SB, the data matrix will often be presented 
graphically by calculating and plotting lines of constant integral values 128 based on 
the individual discrete integral values and their associated positions within die data 
matrix. In some embodiments, this information can be summarized by presenting a 
single line 130 of zero integral value between a region of positive integral values 132 
and a region of negative integral values 134. Alternatively, 3-D or chest anatomy- 
based representations of the data matrix and/or color plots may be used to represent the 
data matrix wherein different colors represent different integral values or ranges of 
integral values. Of course, contour maps or color plots of potential values that 
represent P wave data may also be plotted as set forth above. 
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Referring again to FIG* 2, the P wave data maps are classified visually at step 33 
according to a correspondence in the spatial location and mutual orientation of its map 
extremes and zero line configuration. This will result in one or more groups of P wave 
data maps with nearly identical P wave morphology. Alternative methods to classify 
and group the P wave data maps, eg. pattern recognition algorithms, neural networks, 
statistical analysis or an anatomical subdivision of the location of the left atrial ectopic 
origin, may also be used. At step 34, quantitative data verification is carried out by 
subjecting the groups of P wave data maps to a mathematical analysis using correlation 
coefficients in order to c om p are the uniformity of the P wave pattern within each 
group and the variability of the P wave pattern between groups. Other methods of 
quantitative data evaluation, statistical or the like, may also be used. 

Mean maps are then computed at step 36 to form for instance mean P wave integral 

J\ksc** [ 0 n2ap3 fron a g^up of nearly identical P wave integral maps obtained in different 

ft — • 

\i w pattateJFIG- 6, depicts a group of nearly identical P wave body surface integral 

/StlS^ maps obtained during pacing at the left upper(site A610,B620,D«0,E650andF 
YH^k ; 660) and left lower pulmonary vein (site C 630). The displayed map patterns are all 
ri^oWcife spatially compatible with regard to the location and orientation of the extremes and 
zero hue contour. 

Ja tetyj^i' Referring again to FIQ. 2, the mean P wave data maps are stored in reference database 
VmMt-cft 38. Database 38 contains mean P wave maps { of a variety of ectopic origins in the left 
atrium. The ectopic origin of each mean map is determined for left atrial pacing as 
well as for spontaneously occurring or induced left atrial arrhythmias. In case of left 
atrial pacing, biplane fluoroscopy or other imaging modality is used, such as e.g. 
ultrasound, to determine the anatomical location of the pacing that was performed with 
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a probe or catheter. In case of spontaneously occurring left atrial arrhythmias* similar 
imaging techniques can be used, while multipolar catheter electrogram recordings or 
alternative methods are used to determine the origin of the arrhythmias. 

P wave body surface ECG integral maps of a given patient's left atrial airhythmia 42 
am snbs^iPffitly ^ cnmpared at step 44 to the reference database 38 to classify and 
localize the origin of that patient's left atrial arrhythmia. Consequently, a unique 
database is available for improved classification and localization of left atrial 
arrhythmias The resolution of the database can always be updated and improved by 
addin g new and/or more fax****** information related to the P wave data and the 
underlying ectopic origin sites. 

Referring to FIG* 7, an exemplary database of 17 different mean P wave integral map 
patterns is shown generated by left atrial pacing. Higher or lower number for die 
groups with nearly identical P wave morphology can also be used. A higher number 
of groups improves the resolution of the database and thereby of the classification and 
localization of die left atrial arrhythmias. The encircled numbers 710 (710 A, 710B, 
and 710Q relate maps to a specific endocardial segment of pacing origin shown in the 
anatomical diagrams 720 (720A, 720B and 720C). Endocardial segments of pacing 
can also be delineated as single points, clusters of points, or the Eke. Other anatomical 
or schematic diagrams of the left atrium can also be used Representation of the 
endocardial segments of pacing may also include biplane fluoroscopic views. The 
mean P wave integral map display features extreme positions and zero line contour 
without positive and negative contour lines. Different forms of map format (e.g. 3-D 
or chest anatomy-based format) or map display (e.g. use of various color schemes) can 
also be used The three dimensional diagrams offer a postero-anterior (PA, 720A in 
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FIG. 7A), an anteroposterior (AP, 720B in FIG. 7B), and a left posterior oblique view 
(LPO, 720C in FIG. 7Q of the left atrium. Major anatomical landmarks are 
highlighted: the left (LPA) and right (RPA) pulmonary artery; the superior (SVC) and 
inferior (IVQ vena cava; the left atrial appendage (LAA); right atrium (RA); coronary 
sinus (CS); the left ventricle (LV); the left upper (LUFV), left lower (LLPV), right 
upper (ROPV), and right lower (RLPV) pulmonary vein; Badxmann's bundle (BB); 
mitral annalus (MA); anterior mitral valve leaflet (AMV); and the fossa ovalis (FO). 

Referring to FIG- 8, the design rf an apparatus 140 for dassiiying and localizing left 
atrial arrhythmias is described. The apparatus 140 obtains data for development of a 
database by recording unipolar or bipolar ECO tracings from a plurality of 
predetermined torso sites on a subject 150. Id a preferred embodiment, die tracings are 
obtained from an array of 62 unipolar electrodes during spontaneously occurring or 
induced left atrial arrhythmia or left atrial pacing. Alternatively, a smaller or larger 
number of electrodes or a different distribution of electrodes on the torso can also be 
used 

Electrocardiographic data from the electrodes axe received by a signal processor 152. 
The electrocardiographic data may be acquired simultaneously by the signal processor 
152 at a 1000 Hz sampling rate using a multi-lead ECG recording system. 
Alternatively, the signal processor can perform sequential sampling. Other data rates 
may also be used. The data can be upsampled using multirate filter banks when a 
lower sampling rate is used. Hie signal processor 152 may delete absent signals due to 
electrode obscurement by defibrillator patches or lead dislodgment. The signal 
processor 152 may automatically identify and reject poor quality signals. 
Alternatively, poor quality signals may be visually identified and rejected. These 
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rejected signals may be substituted by interpolating adjacent lead recordings. The 
signal processor may also correct for interelectrode offset differences and linear 
baseline drift, e.g., using linear interpolation after selecting one or more isoelectric 
time instant s or intervals (other techniques for baseline adjustment can also be used). 
The signal processor 152 may incorporate digital filtering to remove line frequency 
interference when necessary, e.g. r with a 60-Hz notch filter (other types of filter or 
means of filtering may also be used). 

A computing means 154 coupled to the signal processor 152 computes one or more 
maps 153, e.g., potential maps, of the P wave, e.g. f at 2-ms intervals (other intervals 
can also be used). The potential maps are visually evaluated to determine P wave 
onset and offset as the point in time when one of the extreme values supersedes and 
returns within a 30 fiV reference voltage window around the zero line. Alternatively, 
other criteria for P wave onset and offset can also be used as well as the onset and 
offset can be automatically detected. In a particular embodiment, an integral map 153 
of die selected time interval in the P wave cycle is computed. As an alternative, the 
computing means 154 may compute P wave potential maps at single or multiple 
selected time instants. The computing means 154 may include software or hardware 
or some c om bi n ation of both. Furthermore, the computing means may be incorporated 
into signal processor 152. 

A validating means 156 then validates the patterns in die maps 153. The validating 
means typically includes circuitry or software for waveform analysis or filtering to 
determine if a given map pattern makes sense. For example, the validating means may 
reject any map fear which the integral or potential values lie outside of a predetermined 
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rang?. Alternatively, the validating means 156 may display the map for visual 
inspection by an expert 

A classifying means 166 classifies one or more groups of P wave data maps with 
nearly identical P wave morphology. Alternative methods to classify and group the P 
wave data maps, e.g. pattern recognition algorithms^ neural networks, statistical 
analysis or an anatomical subdivision of the location of the left atrial ectopic origin, 
may also be used. Once the maps 153 have been classified in groups based on pattern 
resemblance, a verifying means 165 win mathematically evaluate group selection 
using correlation coefficients to compare intragroop pattern uniformity and intergroup 
pattern variability. The verifying means 165 may also perform statistical analysis on 
die map patterns in the groups 153, eg. to determine if the map patterns lie within a 
statistical margin for error. Subsequently, an averaging means 167 averages the 
classified maps that lie in the same classification. The mean maps are then stored in 
the reference database 170. Classified map patterns from a plurality of subjects 150 
build up the database 170. Once the database contains a sufficient number of mean 
map patterns, the apparatus 140 may be used to classify and localize 172 a left atrial 
arrhythmia in a patient 151 by comparing 171 P wave map patterns 153B from the 
patient with the mean P wave map patterns in the database 170. The P wave maps 
153B (e.g. potential or integral maps) of patient 151 are obtained in a similar way as 
for the generation of die database 170 as shown in FIG* 8. This means that 
electrocardiographic data from die patient 151 is reed ved by a signal processor 152B , 
and a computing means 154B to compute one or more maps 153B. 

Referring to FIG* 9, an embodiment of the invention can be implemented as computer 
software in the form of computer readable code executed on a general purpose 
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computer such as compute* 201, or in the fonn of bytecode class files running on such 
a computer. A keyboard 210 and a mouse 211 are coupled to a bi-directional system 
bus 218. The keyboard and mouse are for introducing user input to the computer 
system and c om m u n ic a t in g that user input to processor 213. Other suitable input 
devices may be used in addition to, or in place of, the mouse 211 and keyboard 210. 
I/O (input/output) unit 219 coupled to bi-directional system bus 218 represents such 
I/O elements as a printer, A/V (audio/video), I/O, etc. 

Computer 201 includes a video memory 214, main memory 215 and mass storage 212 
(database storage), all coupled to bi-directional system bus 218 along with keyboard 
210, mouse 211 and processor 213. The mass storage 212 may include both fixed and 
removable media, such as magnetic, optical or magnetic optical storage systems or any 
other available mass storage technology. Bus 218 may contain, for example, thirty- 
two address lines for addressing video memory 214 or main memory 215. The system 
bus 218 also includes, for example, a 32-bit data bus for transferring data between and 
among the components, such as processor 213, main memory 215, video memory 214 
and mass storage 212. Alternatively, multiplex data/address lines may be used instead 
of separate data and address lines. 

In an embodiment of the invention, any suitable microprocessor or microcomputer 
may be utilized as processor 213. Main memory 215 may be a random access memory 
(RAM), dynamic random access memory (DRAM) or an equivalent The main 
memory 215 typically contains a pr o g r am 230 having instructions executable by the 
processor 213. In a preferred embodiment, the program 230 contains processor 
executable instructions that implement die steps of the method 1 of FIG. 1 or method 
21 of FIG. 2. Alternatively, the program 230 may be stored wholly or partly in the 
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mass storage 212. Hie mass storage 212 may contain database 232 of mean P wave 
maps for use with the program 230. The database may have features in common with 
the database 38 and 170 described above in relation to FIG. 2 and FIG. 8, 
respectively. 

Video memory 214 may be a dnal-ported video random access memory. One port of 
the video memory 214 is coupled to video amplifier 216. The video amplifier 216 is 
used to drive the display 217 6^ CRT, LCD). The video amplifier and the monitor 
may be incremented by 

Computer 201 may also include a communication interface 220 coupled to bus 218. 
Communication interface 220 provides a two-way data communication coupling via a 
network link 221 to a local network 222. For example, if communication interface 220 
is an integrated services digital network (ISDN) card or a modem, communication 
interface 220 provides a data communication connection to the corresponding type of 
telephone line, which comprises part of network link 221. If communication interface 
220 is a local area network (LAN) card, communication interface 220 provides a data 
communication connection via network link 221 to a compatible LAN. Wireless links 
are also possible. In any such implementation, communication interface 220 sends and 
receives electrical, electromagnetic or optical signals that carry digital data streams 
representing various types of information. 

Network link 221 typically provides data communication through one or more 
networks to other data devices. For example, network link 221 may provide a 
connection through local network 222 to a local server computer or host 223 or to data 
equipment operated by an Internet Service Provider (ISP) 224. ISP 224 in turn 
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provides data communication services through the .world wide packet data 
communication network now commonly referred to as the "Internet" 225. Local 
network 222 and Internet 225 both use electrical, electromagnetic or optical signals 
that carry digital data streams. The signals through the various networks and the 
signals on network fink 221 and tfarongh communication interface 220, which cany die 
digital data to and from computer 201, are exemplary forms of carrier waves 
transporting die inf ormati on* 

Computer 201 can send messages and receive data, including program code, through 
die netwosk(s) t network link 221, and communication interface 220. In the Internet 
example, remote server computer 226 might transmit a requested code for an 
application program or specific electrocardiography database through Internet 225, ISP 
224, local network 222 and communication interface 220. Consequently, they provide 
an important link for rapid medical information exchange such as telemedicine. 

The received code may be executed by processor 21 3 as it is received, and/or stored in 
mass storage 212, or other non-volatile storage for later execution. In this manner, 
computer 201 may obtain application code in the form of a carrier wave. 

Application code (Le, database, procedures, etc.) may be embodied in any form of 
computer program product. A computer program product comprises a medium 
configured to store or transport computer readable code, or in which computer 
readable code may be embedded. Some examples of computer program products are 
CD-ROM disks, ROM cards, floppy disks, magnetic tapes, computer hard drives, 
servers on a network, and carrier waves. 
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The computer systems described are for purposes of example only. An embodiment of 
the invention may be implemented in any type of computer system or programming or 
processing environment* 

The ability to rapidly acquire, store, analyze and share medical information is 
necessary considering the rapid rate of change in medical technology. This computer 
201 embodiment provides a database system capable of acquiring, analyzing, 
displaying and storing patient data in real time locally* In addition, communicating 
with other medical research centers allows sharing of medical information which can 
add new techniques and enlarge the subject database by the addition of new data 
elements from other centers. 
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